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1 Introduction 

A crucial feature of the behavior of real spin glasses is the existence of extremely 
slow relaxational processes. Any field change causes a very long-lasting relaxation 
of the magnetization and, the response to an ac excitation is noticeably delayed. 
In addition, the characteristics of this slow dynamics evolve during the time 
spent in the spin-glass phase: the systems age. 

Aging effects in real spin glasses have been layed down by experiments 
[1, 2, 3] at a time where there was already an intense theoretical activity on 
the equilibrium properties of mean- field spin-glass models [4] . Experimentalists 
started comprehensive studies of the non-equilibrium dynamics, which happened 
to bring very instructive surprises, while in the meantime theoreticians developed 
extremely sophisticated methods for progressing towards solutions of the equi- 
librium mean-field problem, thence inventing incentive tools for the statistical 
mechanics of disordered systems. 

This early epoch was not the time for the most productive dialogue between 
both parts. The situation is very different now; experiment and theory have had, 
during these last years, a fruitful interplay. On the one hand, the problem of the 
non-equilibrium dynamics has now been theoretically addressed from very dif- 
ferent points of view; scaling theories of domain growth [5, 6, 7, 8], a phase-space 
approach motivated by the Parisi solution to mean- field models [9], a percola- 
tion like picture in phase space [10], random walks in phase space [11, 12, 13], 
mean-field treatments of some simplified situations [14, 15], are now providing us 
with various (and sometimes contradictory) lightings of the experimental results, 
together with impulsing a thrilling debate on the sound nature of the spin-glass 
phase. On the other hand, more and more complex experimental procedures 
[16, 17, 18] have been conceived with the aim of evidencing the materialization 
of some abstract theoretical notions, like e.g. the ultrametric organization of 
states or the chaotic dependence of the spin-spin correlation function on tem- 
perature. 

In this paper, we recall some important experimental features of the spin 
glass dynamics. Since we intend to picture some aspects of the present state of 
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the dialogue between experimentalists and theoreticians, we give a detailed de- 
scription of several ways of scaling the data and of the connection between these 
scalings and the theoretical predictions. We obviously give up any pretention 
of giving an exhaustive coniparison of theory and experiment; wc mainly focus 
here on a perspective of spin glasses which proceeds from mean-field results [14] 
(abundant discussions of the scaling theories can be found in the literature of the 
past few years). In several occasions, we use; as a guideline for the description of 
the experimental results a probabilistic model that views aging as a thermally 
activated random walk in a set of traps with a wide distribution of trapping 
times [12]. 

The slow dynamics of spin glasses - and, as well, of structural glasses [19] 
and other disordered systems - has been often interpreted in terms of thermal 
activation over barriers. One likes to think of a complex free-energy landscape 
(due to frustration) with peaks and valleys of all sizes. This picture has been 
extensively used in the litterature; Refs. [9, 11, 12, 13, 18, 20] are examples of 
different ways of drawing conclusions from it. In fact, the experiments never 
directly probe free-energy valleys or mountains, but rather give access to relax- 
ation rates at various time scales, which may then be interpreted in terms of 
thermal activation over free-energy barriers [18]. 

However, when trying to describe the slow dynamics and aging of real spin 
glasses with a "phase-space" viewpoint, it is worth noticing that phase-space 
is infinite dimensional irrespectively of the finite or infinite dimensionality of 
real space. The geometrical propcirties of the infinite-dimensional phase space 
may put at work a different (non-Arrhenius) mechanism for slow dynamics that 
leads to slowing-down and aging even in the absence of metastable states [21]. 
The particle point in phase space, representing the system c;onfiguration, slowly 
decays through almost flat regions. This mechanism seems to be the one acting 
in the dynamics of mean-field spin-glass models (with a single aging correlation 
scale) [14] as well as in domain growth. In these models one docs not sec neither 
a severe change of behavior when approaching the zero-temperature limit nor 
rapid barrier- crossings from trap to trap in numerical simulations [15, 21]. 

We describe below some mean-field predictions which compare rather well 

with the experiments at constant temperature. Whether one can describe more 
subtle experimental results such as temperature variation dependences, etc. with 
mean-field models and/or with the above non-Arrhenius phase space geometrical 
description is still an open question. The rather good agreement at constant 
temperature suggests that the phase-space dynamical mechanism at the base of 
the dynamics of spin glasses may be a combination of rapid activated processes 
and slow decay through flat regions [14, 21, 22]. We might thence be led to revise 
our "common sense understanding" of the slow dynamics in disordered systems. 
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2 Experimental Evidence for Non-Stationary Dynamics 

2.1 Magnetization Relaxation in Response to a Field Change 

In a measurement of the relaxation of the "thermo-remanent magnetization" 
(TRM), the system is cooled in a small field from above Tg down to some To < 
Tg-, it then "waits" in the field at To during a time tw, after which the field 
is cut, and the subsequent decrease of the TRM from the field-cooled (FC) 
value is recorded as a function of t. Following an "immediate fall-off' of the 
magnetization (depending on the sample and on temperature, of the order of 50 
to 90 %), a slow logarithmic- like relaxation takes place; it is believed to head 
towards zero, although never reaching an end at laboratory time scales. 

These endless-like relaxation processes and, more crucially, the existence of 
"aging" phenomena [1, 2, 3] are a salient feature of spin-glass dynamics: for 
different values of the waiting time t^, different TRM-decay curves are obtained, 
as is evidenced in Fig. l.a. 
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Fig. 1. a. Thermo-remanent magnetization M, normalized by the ficld-coolcd value 
Mfc, vs. t{s) (logiQ scale) for the Ag : Mn2.e% sample, at T = 9K = 0.87Tg. The 
sample has been cooled in a 0.1 Oe field from above Tg = IQAK to 9K; after waiting 
tju, the field has been cut at t = 0, and the decaying magnetization recorded. 



The dynamics depends on two independent time-scales, t ( "observation time" ) 
and tw ("waiting time"). This dynamics is non- stationary: the response at t + t^ 
to an excitation at tyj depends on t + ty, and t^, and not only on t (breakdown 
of time-translational invariance). Qualitatively, one can see in Fig. l.a that the 
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longer the waiting time before cutting the field, the slower the overall response; 
the initial fall-off is smaller, the relaxation curve shows a slower decrease, the 
system has become "stiffer". Such aging phenomena have been early identified 

in the mechanical properties of glassy polymers ^ [23]; the slow strain following 
the application of a stress has been recognized to depend on the time spent in 
the glassy phase. 

In spin glasses, the aging phenomena have initially been explored using the 
mirror experimental procedure of the TRM [1], in which the sample is cooled 
in zero-field, and after tw a small field is applied. As far as the field remains 
low enough (usually, in the range 0.1-10 Oc), both procedures are equivalent; 
the relaxation of this "zero-field cooled magnetization" (ZFC) follows the same 
t and dependence as the relaxation of the TRM [24, 25]. More precisely, 
it has been shown in [24] that, for all t^tw values (i.e. all along the measured 
relaxations for various tw), the sum of the ZFC-magnetization plus the TRM 
equals the field-cooled value. This is simply linearity in the response, since this 
experimental result shows that the sum of the responses to different excitations 
is equal to the response to the sum of both excitations (the response to a constant 
field being the field-cooled magnetization). That linearity holds for all t,tw tells 
us that the presence of the (sufficiently small) field does not influence the aging 
process: waiting t^ in zero field and then applying a field during t (ZFC case) 
is equivalent, for the dynamics, to applying a field during tw and then waiting t 
in zero field (TRM case). The only role played by the field in this context is to 
reveal the dynamic properties of the system. A recent study of the effect on the 
dynamics of increasing field values can be found in [20, 26, 27]. 

In the semi-log plot of Fig. l.a, each curve shows an infiection point, and 
one first quantitative estimate of the -effect on the relaxation is that this 
inflection point is located around log t ~ log tw ■ This fact has been noticed and 
given a physic;al meaning by Lundgren et al. [1]. The relaxations are slower than 
exponential; they do not correspond to a single characteristic response time r, 
but are likely to be parametrized with the help of a wide distribution gt^ir), 
which is defined hereby: 



where tq ~ 10~^^sec is a microscopic attempt time. M{t + tw,tw) is the ZFC or 
TRM, depending on the experiment, and Mfc{t + tw) is the field-cooled value at 
time t + tw In the figures we abbreviate M(t + tw,tw)/Mfc{t + t^) = M/Mfc- 
Lundgren et al. have pointed out that taking the derivative of (1) with respect 
to log* gives access to the distribution gt„{T), since 



^ A wide class of materials like e.g. PVC, PS, Epoxy, or even bitumen, Wood's metal, 
amorphous sugar and cheese [23]. 




(1) 
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a rough approximation reflecting the sharp character of ^ exp(— |:) around t = t 
(to be considered on a logarithmic scale, which is actually the scale which is 
suggested by the measurements). The plot of the relaxation derivatives shows 
bell-like shapes, with a broad maximum around logt = logtyj, and pictures 
gt^ir) for various tw [1]. Thus, in a first approximation, the aging phenomenon 
can be described as a logarithmic shift towards longer times of a wide spectrum 
of response times ^. This shift is of the order of log tyj, and therefore suggests 
that the dynamics be the same as a function of t/t^. 

Let us call "full aging" the pure t/tw scaling, that is not far from being the 
correct one, as seen in Fig. l.b where the data from Fig. l.a is presented versus 
t/tyj. Most of the t^-effect has been accounted for, though some systematic 
departures remain, and are worth being discussed (see Sect. 4.2). 




2.2 Ac susceptibility 

The approximate t/tu, scaling of the TRM (or ZFC) curves is sufficient for a 

description of aging effects in ac experiments, where the in-phase and out-of- 
phase components of the response to a small ac excitation field at a frequency 

^ Indeed, the spin-glciss properties do not exactly depend on tw, but rather on the total 
elapsed time tw+t (Sect. 4.2); they are evolving during the TRM measurement itself 
[3, 28]. The physical interpretation [1] of gt^ir) therefore remains approximate. 
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to are measured. Aging is more visible (in relative value) in the out-of-phase 
component x" of the magnetic susceptibility, which represents dissipation. The 
observation time, corresponding to t in TRM experiments, is here constant, 
equal to l/uj. When the sample is cooled from above Tg down to Tq < Tg, 
the susceptibility does not immediately reach an equilibrium value, but shows 
a slow relaxation as time goes on. We denote ta ( "age" ) this time elapsed from 
the quench into the spin-glass phase; in the TRM experiment, the equivalent age 
is t + tyj = ta- The non-stationary character of the dynamics, which appears in 
the TRM measurements as a dependence on the two independent time scales t 
and tw, shows up in ac experiments as a dependence of x" on the two variables 
w and ta- This is clear in Fig. 2 where x" at various (low) frequencies iv is plotted 
as a function of cot a; applying a vertical shift, the curves can all be merged with 
respect to this reduced variable, which is equivalent to t/t^ in TRM experiments. 
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Fig. 2. Out of phase susceptibility xi^j^a) vs. ujta for the insulating CdCri.rlno.sSi 
sample. The four curves, corresponding to lj = 0.01,0.03,0.1, 1. Hz, have been verti- 
cally shifted (see text), ta is the total time elapsed from the quench (age). The inset 
shows a scaling of the same data which follows from mean-field results (see Sect. 3.2). 

Thus, the approximate t/tw scaling obtained from TRM and ZFC experi- 
ments can be fairly well transposed to an co.ta scaling of the ac susceptibility 
X"{ui.ta). The vertical shift corresponds to accounting for the various "equilib- 
rium values" (x"eq(w) = limx"(w,ta ooj) at different frequencies. In Fig. 2, 
for technical reasons, the zero of the scale has not been measured and the shift 
is arbitrary. As an example of the relative orders of magnitude, let us mention 
that, at Lo = O.OlHz, the amount of the relaxing part is roughly equal to the 
equilibrium value; for %' in the same conditions, it would be of the order of 10% 
of the equilibrium value. The frequency-dependence of x" eq has been determined 
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in other studies [25, 28]; it can be represented by a power law with a very smaU 
exponent (or else a power law of a logarithm) 

X"eg(t^) OC , (3) 

where a increases in the range 0.01 — 0.1 when approaching Tg from below. This 
is valid in the 10^^ — \{y^Hz range which has been explored, and has been mea- 
sured rather in insulating than in intermetallic spin glasses (due to eddy currents 
in metals). Both classes of samples have been found to present the same gen- 
eral spin-glass behavior [3, 28]. The LO.ta scaling indicates that the smaller the 
frequency, the longer the time ta during which a significant relaxation, character- 
istic of aging effects, can be found. Therefore, at higher frequencies (w > lOifz) 
aging disappears very rapidly, yielding almost instantaneously a stable value 
X"eg(w); conversely, at lower frequencies, the determination of x"eg(w) becomes 
problematic, implying measurements over tens of hours or days. 

2.3 Time Regimes in ac and dc Experiments 

Let us summarize the conditions for observing either stationary or non-stationary 

(aging) dynamics in ac(x") and dc (TRM or ZFC) measurements. We can define 
two distinct time regimes, which apply to both experiments. 

— For Lo.ta ^ "equilibrium dynamics" is recovered in ac experiments, in the 
sense that only one time scale is needed: the dynamics is then stationary. 
This time regime corresponds, for TRM's, to {t^^ + t)/t = ta/t » 1 or 
equivalently t <^tw, that is the very beginning of the TRM-decay curves. 

— If LJta is comparable to ta [iota = 0(1)) in the x" experiments and, equiv- 
alently, t is comparable to t^ in the TRM experiments, one observes non- 
stationary dynamics. 

The x" measurements are limited to uta > 1 since uj.ta < 1 cannot be ex- 
perimentally realized (the harmonic response is not defined at times shorter 
than one period). Hence, with x" we can only explore the beginning of the 
aging regime. 

In contrast, aging is predominant in TRM-measurements over the largest 
part of the accessible time scale since, for TRM's, (cj.ta) = ta/t = {tw +t)/t 
becomes rapidly close to 1 as the observation time t elapses. Measuring the 
TRM decay we have access to a larger time-window in the aging regime. 
Therefore, aging in x" can only be compared with the TRM decay at the 
beginning of the aging regime - that we call early epochs (also called "quasi- 
stationary regime" in [28]). This might explain why a full aging scaling seems 
to apply better to x" results than to the TRM (see Fig. l.b and 2). 
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3 Aging Theories for Old Results 
3.1 Scaling theories 

The different time scales in which stationary and non- stationary dynamics are 
occurring arc likely to be mapped onto length scales in the real space of spins. 
When approaching Tg from above, the onset of a critical regime has been char- 
acterized by the diverging behavior of the characteristic time in ac experiments 
and of the non- linear susceptibility in dc studies, [29]. A thermodynamic phase 
transition would imply the divergence of a characteristic correlation length ^ 
when T ^ Tg. In this equilibrium picture the spin-glass phase is believed to be 
an ensemble of randomly oriented spins, which are frozen due to infinite-range 
correlations corresponding to a long-distance "order" . 

However, aging shows that indeed equilibrium has not been established when 
crossing Tg. In several "scaling theories" [5, 6, 7, 8] of non-equilibrium phenom- 
ena in spin glasses, the spin correlations are considered to be limited to some 
finite range ^{t) (out-of-equilibrium situation); as time elapses, spin rearrange- 
ments yield a slow (due to frustration) extension of equilibrium correlations, 
towards the equilibrium situation of infinite range {t oo,^{t) oo). An ac 
experiment at frequency w, as well as a TRM or ZFC relaxation at time t, can 
be viewed as probing the spin-glass excitations at a given length scale L which 
should be an increasing function of the characteristic probe time 1 /oj for x" and 
t for TRM. In the "droplet model" by Fisher and Huse [6], one has L oc log^^''' t 
{tp < d ~ 1), and in the "domain model" by Kopcr and Hilhorst [7] L oc t^^'^ 
{p ~ 0.5). These theories do not aim at a microscopic description at the scale of 
spins, but make this mapping of time onto length scales quantitative, in terms 
of scaling laws. At least qualitatively, these models provide us with a convenient 
picture of aging phenomena, which is the following. For short probe times com- 
pared to the age {ui.ta ^ 1, or t/t^ <C 1), short-ranged excitations are involved 
(L ^ C(i)), and the increase with time of ^(t) does not affect the dynamics, 
which is found to be stationary (no aging). Conversely, for longer probe times 
compared to the age {u).ta ~ 1, or t > tyj), the characteristic length of the rel- 
evant excitations is of the same order of magnitude as $,{t), which is increasing 
due to aging, and the dynamic properties are strongly affected by aging (non- 
stationary dynamics). 

The accurate quantitative agreement with the data is still to be discussed 
(see e.g. Ref. [13]). Some critical remarks to the simple scaling approaches, which 
rise up in view of other results (T- variation experiments), are discussed in Sect. 



3.2 Non-Equilibrium Dynamics in Microscopic Theories 

The Models. The classical "realistic" microscopic model of spin glasses is the 
3-D Edwards-Anderson (3DEA) model [30] 



5.1. 




(4) 



{iJ) 
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where are Gaussian or bimodal random variables, Si are Ising spins and {i,j) 
represents a sum over first neighbours on a cubic 3D lattice. It is very difficult to 
obtain analytical results for the statics or the dynamics of 3DEA, in consequence, 
even after more than 20 years of research on the field of spin glasses, very few 
results are available. A lot of efforts have been devoted to the numerical study 
mainly of the equilibrium properties of the 3DEA. Again, the situation is still 
pretty unclear: basic questions as to the existence of a thermodynamic phase 
transition are still not answered [31]. 

The standard mean-field extension of the 3DEA model is due to Sherrington 
and Kirkpatrick (SK) [32] and corresponds to the same interactions as in (4) but 
with the sum extended to hold over all pairs of spins in the system. The study 
of the SK model - and of some other related mean-field models - had for a long 
time been confined to the search for equilibrium properties [4]. 

Numerical Simulations of Out-of-Equilibrum Phenomena. One may 

wonder whether aging, which did at first sight appear as some imperfection 
of the experiments, is really intrinsic to the Hamiltonian (4). We now know that 
the answer is yes. Only recently, attention has been paid to the study of the out 
of equilibrium dynamics of microscopic spin- glass models. Andersson et al. [33] 
and Riegcr [34] reproduced in a numerical simulation the procedure of, e.g., the 
TRM experience using the 3DEA model. The results show that it captures the 
main features of real spin glasses: both slow dynamics and aging effects. Later, 
numerical simulations of the large D "hypcrcubic" spin-glass cell in real space 
showed that also this model, that is expected to reproduce the SK model, when 
— > 00, captures the main characteristics of aging [35], thus confirming the 
previous analytical results that we describe in the following paragraphs. 

Analytical Approach: Formulation and Definitions. Again, it only hap- 
pened recently that analytical developments evidenced aging effects in mean-field 

spin-glass models [14], showing that these simplified models can describe, at 
least qualitatively, the phenomenology of real spin glasses. The idea in the case 
of mean-field spin-glass models is just to try to solve the exact dynamical equa- 
tions derived for N, the number of dynamical variables in the model, tending to 
infinity. These equations are well-defined, have a unique solution and, in the ab- 
sence of a magnetic field, only involve the correlation and the response functions 
(sec Eqs.(5),(10) below for their definitions). The initial condition is chosen to 
be random so as to mimic the initial configuration just after the quench in the 
experimental situation. One then considers the large-time limits, but only after 
having already taken the thermodynamic limit A'' — > oo to obtain the asymptotic 
behavior of the solution. 

It is important to notice that in this approach it is not necessary to assume 
a priori any particular structure of phase space - typically to say that there are 
many metastable states due to frustration separated by high barriers - to obtain 
the dynamical behavior of the problem. The solution can a posteriori be given 
a geometrical interpretation [12, 13, 21, 22]. 
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The solution shows that the equations are self-consistently solved in the large- 
time limit by an aging solution. The reason why these equations can be solved 
is the weak long-term memory of the system [14]. The results fall into the weak- 
ergodicity breaking scenario previously proposed in [12, 13] within the trap model 
(see Sect. 3.3 below). When looking at the large-time dynamics of the system, 
the weak long-term memory property allows us to neglect the contribution of 
any finite time-interval after the quenching time. The system forgets what hap- 
pens in "finite" time intervals with respect to the "infinite" observation time. It 
keeps, however, an averaged memory of its history. The weak-ergodicity breaking 
scenario tells us that the evolution of the system continues forever; the dynamics 
slows down as time elapses but the system is never completely stopped in its 
evolution. The waiting-time tw gives us an idea of the age of the system. 

In the following we shall be a bit more technical and describe the main 
features of the formalism and the solution. 

The auto-correlation function is defined as 

1 ^ 

C{t-\-t^,t^) = — 'Y^{si{t + t^)si{t^)) , (5) 

1=1 

with the overline representing a mean over different realizations of the disorder 
and ( ) an average over different realizations of the thermal noise. We then define 

CF{t)= lim Cit + ty„ty,) Cf{0) = 1 \im CFit)=qEA ■ (6) 

tu! — *00 t — >00 

This allows us to separate the auto-correlation into two additive parts, a station- 
ary term and an aging term Ca [14, 13]: 

C{t + t^,t^)=CF{t)-qEA + CA{t + t^,t^) . (7) 

In the absence of a magnetic field, the weak ergodicity breaking scenario [12, 13, 
14] implies 

lim C{t -\- tu,, tu,) = V fixed t^ , (8) 

thus 

hmt^oo lim(„,^oo CA{t + tw,tyj) = qea /q\ 
limt^^CA{t + t^u,tw) = . ^' 

It will turn out that the two scales corresponding to these two limits are well- 
separated for mean-field models, in the sense that in the time-regime where Cf 
varies then Ca stays constant, and viceversa. In other words, one can think of qea 
as a value of the correlation separating different "correlation-scales" , C > qea 
and C < Qea' when t <C t^, C > qEA and we have stationary dynamics, while 
when t^ tw, C < qea and we have non-stationary dynamics and aging just as 
described in Sect. 2.3 for the general features of aging in spin glasses. ^ 

^ In some numerical works, the form C{t -\- tvjjtyj) = t~'^^^^ ${t/tw) has been often 
used to scale the data for all times t ([34, 35, 36, 37], see also [28] for a related 
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In the same way, the response function can be equivalently separated into a 
stationary and a non-stationary term 



= RF{t)+RA{t + ty„ty,) . (10) 

with 

RF{t)= lini R{t-\-tyj,tw) and (11) 

tw—*oo C{t+tw ,tw)>qEA 

RA{t + ty„t^)= lim R{t + t^,t^) . (12) 

Rpit) satisfies the fluctuation-dissipation theorem (FDT) and RA{t + tw,tw) 
satisfies a generalized FDT [14] 

1 dCpit) X[CA{t + t^,t^)] dCA{t + t^,t') ^ 
iiF\t} — - — — KA[t+tw,tw) - ^ \t'=tw 

(13) 

with < X[Ca] < 1 a monotonicahy increasing function of < Ca < Qea- 
{X = 1 corresponds to the usual FDT.) 

With these definitions one can solve the asymptotic (large tw) dynamics of 
several mean-field disordered models [14, 15, 38, 39, 57]. 



1 ^ 

R{t + tw, tw) = ^ 



d{s,{t + t^)) 



^hi{tw) 



Analytical Approach: Stationary Regime. For all these models, when t is 
large, the stationary part of the correlation function Cpit) decays with a power 
law 

CF{t)r^qEA + Ca('jY , (14) 

where tq is a microscopic time-scale and a has precisely the same meaning as 
the exponent in (3), describing the frequency dependence of the equilibrium out 
of phase susceptibility. 

The temperature dependence of a depends on the model. For models ^ such 

as the p-spin spherical spin glass [40, 39] or the model of a particle moving 
in an infinite-dimensional random potential [41, 39], a = 1/2 at T = and 

discussion of the experimental data). It should be remarked that, though at first 
glance this scaling seems to be similar to the one following from the WEB scenario, 
it implies quite a different conclusion for the global behavior of the system. Note 
that if one takes the limit limt^oo limt^^oo, that corresponds to exploring the end 
of the stationary dynamics, Eqs.(6) yield limt^oo limt^^oo C(i-f i™, tu,) = qea, while 
limt^oolimt^^oot~'^^'^^ 4>{t/tw) = *(0) limt^oo i"""^^' = 0. The stationary dynam- 
ics and the fact that the correlation decays to qea at the end of this time-regime, 
have a very clear geometrical interpretation [13, 14, 21]. 
* Though the aim of Refs. [40, 41] was to study the equilibrium dynamics a la Som- 
polinsky [42] - a different situation from the out of equilibrium occuring in exper- 
iments - the calculation of the exponent a obtained in these works applies to the 
experimental case when adequately reinterpreted [14]. Let us also note that the in 
this paper the a and f3 exponents are exchanged with respect to Refs. [39]. 
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it decreases when increasing the temperature. For the SK model, conversely, 
a=l/2atr = Tc and it decreases when decreasing the temperature [43]. 
Finally, for the mixed {p = 2 + 4) spherical model introduced in Ref . [44] , a has 
a non-monotonic dependence on T; a{T = 0) = a{Tc) = 1/2. The value of the 
exponent a measured experimentally follows the tendency of the one holding for 
SK and the mixed (p = 2 + 4) spherical models close to the critical temperature, 
though the value of a from the experiments is considerably smaller (a < 0.1 vs 
a ~ 0.5). 

Analytical Approach: Non-Stationary Regime. Following then very gen- 
eral requirements, it has been argued in Ref. [14. 15], and explicitcly chocked 
on several pure and disordered models, that in the large-time limit only two 
situations with different dynamical behavior seem to exist: 

— On the one hand, there are models with only one time-scale - or equivalently, 
correlation-scale - apart from the stationary one. Ca scales as in a domain 
growth process within the non-stationary time-scale, in the sense that: 



with h{t) a monotonically increasing fimction - analogous to the domain 
length L{t). j~^{u) is a function characterized by another exponent which 
we call here (1 — x) (and has been called /? in [39] and a in [45]). Close to 
u = l, i.e. for the early epochs of the aging regime, reads 



c is a constant and x < 1 implying that is non-analytical in the neigh- 
bourhood of u = 1 (see (35) below). 

In this case, the FDT-violating factor X[Ca] is a constant X < 1. These 
models, when treated statically with the replica trick, are solved by a one 
step replica symmetry breaking ansatz [4] . An example is the p-spin spherical 
model [46]. We call them "single-scale models". 

Certainly the functions h and do depend on the specific model. At the 
mean- field level we have succeeded in obtaining and X for several mod- 
els. However, surprisingly, there are for the moment no analytical results 
available for the scaling function h{t) ^. The simplest possibility is that h{t) 
is a pure power-law (t/ro)", as found in the trap model [12] or standard 
coarsening models. This solution is particular in the sense that (15) is then 
independent of the microscopic time scale tq, which can be taken to zero; in 
this case Ca simply depends on the ratio t/tw (full aging situation, a quali- 
tative approximation of the experimental results, as explained in Sect. 2.1). 

^ This technical difficulty is related to the introduction of a time re-parametrization 
invariance when studying the exact mean-field equations for large and widely sepa^- 
rated times tw and t + tw. 




(15) 



J ^{u) ~ qea - c(l - u) 



1-x 



(16) 
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This is not the case for more general functional forms. Two explicit choices 
which have been proposed so far are: 



The form on the left was proposed to account for experiments in polymer 
glasses by Struik [23] , then used in the first accurate analyses of aging effects 
in the TRM-decay [3, 28], and recently found in the exact solution of the 
asymmetric spherical SK model with = 1/2 [47] (see also [38]). The second 
form is suggested by the numerical data from the "toy model" of a point 
particle in a random potential with infinite dimension [39] . Both will be used 
below to scale the data for the TRM and the out-of-phasc susceptibility. Note 
that in the limit ^ = 1 or a = 1, one recovers full aging with a pure power- law 
behavior for h, while n = corresponds to time translation invariance (no 
aging). When ^ < 1 and a > 1 we have "sub-aging" that we define as follows. 
Taking t fixed and, say, in the beginning of the aging regime, h{tw)/h{t + 
tw) ~ 1 — {d\n{h{tu])) /dtw)t. This defines a characteristic relaxation time 
T{tw)- We say that we have sub-aging (super-aging) when T{tw) grows slower 
(faster) than tw 

Interestingly enough, one can in general derive a relation between a, (1 — x) 
and X (see (14), (16) and (13) for their definitions) [39, 45]: 



for single scale models. We shall use this equation to predict X at the be- 
ginning of the beginning of the aging regime in Sect. 4.2. 

— On the other hand, there are models such as SK that have an infinite number 

of time-scales - correlation scales - apart form the stationary one [42, 15, 38]; 
mathematically, one has ultrametricity in time for all correlations such that 
Ca < qsA, in the sense that (^^(^1,^3) = min(CA(ti, i2), CA(t2, is)), h > 
t2 > ts [15]. The decay is here infinitely slower than in the single-scale 
models. The FDT-violating factor X[C>i] is a nontrivial function of Ca- 
These models are solved by a full replica- symmetry breaking ansatz when 
using the replica trick at the static level [4] , and can be called "multi-scale" 
models. 

It is important to notice that a scaling like (15) inside a correlation scale 
and ultrametricity between different correlation scales are expected to hold on 
very general grounds, in particular for more realistic finite dimensional models 
(in the limit of large-times) , provided that the rather mild assumptions used in 
[14, 15] are satisfied. This justifies the fact that we shall use, in the following, a 
scaling-law like (15) to scale the data for real spin glasses, without refering to 
any particular model. 




or 




{r[i + {i-x)]r _ {r[i-a]r 

r[l + 2(l-a;)] r[l-2a] ' 



(18) 
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Connection with Measurable Quantities. If linear response theory holds, 
the TRM is just 

M{t + t^,t^) = h [ dsRit + t^,s) . (19) 

For large waiting-time tw, this integral can be rewritten using the decomposition 
in stationary and non-stationary decays ((10) and (13)) and using (14), for t ^ tq 
we have 



^(t)"^/ ^^^aXICa]. (20) 



Mfdt + tw) 

where A is a constant. For single scale models as in (15) the scaling reads 



Mfc{t + tw) \tJ \h{t + tn^) 

In the early epochs of the aging regime for the TRM, that should be compared 
to the non-stationary behavior of the x"('^ji)' one has 

Mfcit + tw) 
where we used (16). 



The out-of-phase susceptibility can also be simply related to the correlation 
function (5). For high-frequencies, ut oo, the aging term does not contribute 
(it is the integral of a slowly varying function Ra (t, s) times a rapidly oscillat- 
ing function). Using (14) to approximate the remaining integral, one finds the 
stationary part of the a.c. susceptibility: 

X"{UJ, t) ^ Xe,(^) « , ^ 00 . (23) 

Conversely, \i ujt > 1, i.e. for low frequencies, the aging part strongly con- 
tributes. For single-scale models we then have 

X"M - XU^) ~ ^huj^ dseMio^s) , u,t>l. (24) 

For uit finite but large, one has in general: 

x>.*)-.»«r'(i)-r'(i-^^) (25) 

fd\nh(t) iV'"' 
"^[-^Z^t) l«-*<-, (26) 

where we introduced the power- law behavior of J~^{u) in the vicinity of u = 1 
defined in (16). Hence the conclusions for x": 
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— If h{t) is a simple power-law, then x"{^j t) ~ x'egi^) scales as cot, as obtained 
in [12, 13] (full aging). 

— If h{t) = exp(l/(l - m)*^"^) then x"(w,t) - x",(w) is a function of cot x 
{t/To}^^^- When n 1 there is a correction to the pure uit scaling (sub- 
aging for fj, < 1)). 

— If h{t) is of the form exp[ln°(t/To)], then the wt scaling is corrected by a 

slowly varying factor In'-'^^^-' {t/ro) (sub-aging if a > 1). The x" data are 
scaled within this assumption in the inset of Fig. 2. 

In Sect. 4, we apply the scaling relation (15) from single-scale models together 
with these proposals for h{t) to the TRM and x" data. 

The mean-field models which lead to the above results are very instructive: 
general statements and new ideas (like the violation of FDT) have emerged from 
their study. However, the physical mechanism underlying aging in these models 
is not yet very clear. The single-scale models, in particular, are very weakly 
sensitive to temperature, suggesting that no activated effects are involved, and 
that aging is rather related to large dimensional effects: the system wanders 
indefinitely in a large phase-space, without ever reaching a local minimum of the 
free energy [21]. Conversely, the trap model which we shall discuss now relies on 
activated effects to generate a broad distribution of time scales, which also leads 
to aging. However, this model is phenomenological, and does not emerge from 
a precise microscopic description - although some steps in this direction have 
recently been made [22, 48]. A tentative classification of the different models of 
aging has been proposed in [49]. 

3.3 The Trap Model 

In the trap model [12], aging has been shown to naturally occur in a situation 
called "weak ergodicity breaking" , which corresponds here to a statistical impos- 
sibility for the system to realize equilibrium occupation rates of the metastable 
states. This model has appeared as a fertile guideline for the analysis of the ex- 
periments; we therefore recall its main points, and come back to it later in Sect. 
5.2. In the simplest version of the model [12], aging is sketched by a random walk 
in a collection of "traps" with random trapping times t, all equally accessible. To 
each trap is associated a certain magnetization M and ac susceptibility Xt{i^)- 
The properties of a real sample arc obtained by averaging over an ensemble of 
decorrelated subsystems, corresponding to spins in different regions of space ^. 
An important input is common to microscopic theories [4] , and also to the more 
general problem of manifolds in random media [22, 48]; the distribution of trap 

^ At this stage, the spins do not enter directly. In later developments [26], however, 
the number of spins to be flipped for escaping from a trap (and hence the size of 
the subsystems referred to above) has been estimated from the influence of the fleld 
amplitude on the dynamics, as observed in the experiments. 
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depths is taken as an exponential. For thermally activated processes, this yields 
the following distribution of trapping times: 

V'W = :^ (forr»To) , (27) 

where x (from the distribution of barrier heights) is a temperature dependent 
parameter describing the structure of the phase space. In a comparable way, the 
"random energy model" (REM) of Dcrrida [50] involves x = T/Tg. The crucial 
point is that a; < 1 in the spin-glass phase; in consequence, the mean value of ^/'(t) 
is divergent, that is the mean time needed to explore the whole set of traps (and 
thus to reach crgodicity) is infinite. This was called weak crgodicity breaking, in 
the sense that the equilibrium situation is never realized, although the system 
never gets trapped in a finite region, leading to an asymptotically zero correlation 
function - sec (8). This is very different of the usual crgodicity breaking, where 
the system can reach rather quickly an equilibrium configuration, but remains in 
a restricted sector of the phase space. Since x <1, the distribution (27) is very 
broad. After a random walk during the system has visited numerous short-life 
traps, but in ci rclcLtivcly sina,ll time coinpcirccl with. tyj\ as usual with such broad 
distributions, the significant contributions arise from the largest - although rare 
- events. Thus, after t^, the system has the largest probability to be found in 
a trap of characteristic time itself; if a magnetic field is varied at t^, most 
subsystems will need a time of order before changing their magnetization. It 
has been shown [12] that the TRM-decay is then a function of t/t^, and similarly 
that the ac susceptibility is a function of uj.t. Thus, on the basis of a statistical 
description in the space of the metastable states, the main features of aging can 
be obtained. Further developments of this approach [13, 52] are discussed below 
at the light of various aspects of the experimental results. We now turn to a 
more detailed description of the combination of aging and stationary dynamics 
as seen in both TRM and x" measurements. 



4 TRM Experiments: Aging and Non-Aging Dynamics 
Disentangled 

4.1 Departures from Full t/tw Aging 

In Fig.l.b, the TRM curves are presented as a function of t/tw] it is clear that 
this is not exactly the correct reduced variable (failure of a full aging scaling). 
The same effect has been found in other samples [3, 28] and can be seen in 
results from other laboratories {e.g. [1]). It has been first identified in polymer 
mechanics [23]. Indeed, the x" results recall us that stationary dynamics, namely 
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X"eg('^) oc (an additive contribution to the aging part, see Fig. 2), must 
intervene in the TRM decay, particularly in the t <^tw regime. This frequency- 
dependent x"eq(t^) is equivalent to an additive contribution of the form t~°' to 
the TRM. 

The question of the departure from full aging has already been discussed in 
the past and the TRM's have been very accurately parametrized as the product 
of a factor times a {t,tw) dependent factor ([28], see also the footnote^ 
above). However, an additive combination of the stationary and aging parts arises 
naturally in the above theoretical approaches [14, 12], and we shall reanalyze 
these data in this light. We express the stationary part of the TRM, in units of 
the field-cooled value Mfc like in Fig. 1, as A{To/t)°'. tq is again a microscopic 
time, which allows homogeneity of the units, and then A is a non-dimensional 
constant, expected to be of order 1. We have taken the same TRM-data as in 
Fig. 1, and adjusted A and a in order to try to merge the aging parts f{tyj, t) 




of all 5 curves of various t^s as a function oit/tw We have no a values from x"- 
measurements on the AgMn sample, so we have simply kept a in the 0.01 — 0.1 
range obtained for the CdCri^^Ino.sSi sample in previous analyses [28]. Fig. 3. a 
shows the result; whatever the choice of parameters (also, trying a logarithmic 
decrease instead of a power law), it seems impossible to merge all 5 curves 
together in the whole time regime. Systematic discrepancies are always found 
in the large t/tw range: no simple function of t alone can account for these 
deviations. However, for small t/tw, the good quality of the scaling is equivalent 
to the simple w.t-dependence of x" (w, t), which is actually measured in the same 
time regime (i.e. cot > 1). 

In the long time t/t^ > 1 region, the discrepancies shown in Fig. 3. a present 
the same features - although somewhat weaker - as in previous analyses using 
a multiplicative stationary contribution [28]. In Fig. 3. a, as well as in Fig. l.b 
(the same, without stationary dynamics subtraction), it is clear that a t/tw 
scaling slightly overestimates the aging effect (sub-aging); the "youngest" curve 
(shortest tw) lies above the others, whereas it was below in the raw data (Fig. 
l.a), and the effect is sytematic for the 5 curves. 

4.2 Sub-Aging Scaling 

A way to succeed in merging all TRM curves is to use a generalized scaling 
function of the form h{t) = exp[Y^(;^)^^^] with ji < 1 (one of both examples 
quoted above in (17)). This was proposed in the context of spin glasses in [3, 
28], as a phenomenological scaling procedure inspired from polymer mechanics 
[23], which proved to account with great accuracy for the observed sub-aging 
situation. We refer the reader to [3, 28] for the arguments which lead to the 
above choice of h{t), or rather, along the lines of [28], to the effective time A 
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Fig. 3. a. Aging part of the TRM (Eq.28): the estimated stationary contribution has 
been subtracted from the frill measured value. The data (same as in Fig. lab) is plotted 

vs. t/tw 



defined as 



7 . 



As a function of A, the 5 curves in Fig. 3. a recover precisely the same shape 
[3, 28], as displayed in Fig. 3.b. In other words, A results from a change of 
variable which allows us to see the aging dynamics as stationary. It can also be 
obtained from 

which means that, since the age tyj +t is varying during the TRM experiment, 

the effect of an elementary time interval dt depends on the elapsed time; due to 
sub-aging, the effective time scale associated to the age (tw+t) is {tyj + t)'^, with 
At < 1. 

This type of scaling has been successfully applied to various samples [3, 28]; n 
is found in the 0.8 — 0.9 range and it is almost independent of temperature in the 
0.3 < T/Tg < 0.9 range. This /Li-trick is a very convenient way of parametrizing 
the sub-aging deviations from a full i/t^-scaling, although not necessarily having 
a direct physical meaning, /x = corresponds to the case of no fu,-dependence 
(no aging), while /U = 1 would yield a full t/tw scaling. The intermediate values 



^ The careful reader will notice a change of A-units when compared with [28]; within 
the present definition, A is equal to from [28]. 
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Fig. 3. b. Aging part of the TRM as in Fig. 3a, but vs the scaling variable A/t^ defined 
in (29) (1st example in (17)). 5 curves, obtained for tyj = 300, 1000, 3000, 10000, 30000s 
(same data as in Fig. 1 and 3a), are superimposed onto each other. 



of n (sub-aging) correspond to the fact that the apparent relaxation time scales 
sub-linearly (as {tyj + tY) with the age + t. 

The alternative choice of h{t) proposed above in (17) can also be considered 
[53]. It is actually very close to the /i-case, in particular in the limit /i — > 1. 
Again, for a = 1, full aging is recovered, while a > 1 describes a sub- aging 
situation, with an apparent relaxation time T{tw) = tyj/{a\n^°'~^\t/To)). We 
have also applied this other sub-aging scaling to the non-stationary part of the 
relaxation displayed in Fig. 3a; the result is shown in Fig. 3.c, again the curves 
are fairly well superimposed onto each other. 

Eq.(18), derived for single-scale mean-field spin-glass models, relates a, x 
and the FDT violating factor X defined in (13). According to (15) and (16), we 
have made a fit of the beginning of the TRM decay in Fig. 3.c (early epochs), 
and we have obtained x ^ 0.96 (sec (16)). This is consistent with X = 1, which 
would suggest that FDT applies without corrections even at the beginning of the 
aging regime, if one accepts that (18) holds. This would be in accord with previ- 
ous indications obtained through the comparison of x" and noise measurements 
made in the early epochs of the aging regime ( "quasi-stationary regime" ) [28] . A 
careful and detailed study of the noise auto-correlation and response functions 
would help us knowing if and how FDT is violated, namely if X departs from 1 
when going deeper in the aging regime. . 

For the sake of comparison, we show in Fig. 3d how the /x-scaling applies 
to the full value of the TRM (the stationary contribution is neglected, and not 
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Fig. 3. c. Same as Fig. 3b, but as a function of a scaling variable which corresponds 
to the 2nd example in (17). 



subtracted from the magnetization). An almost acceptable scaling is obtained; 
the discrepancies are not much larger than the experimental uncertainties, but 
they are systematic: all curves are crossing in the middle of the figure. This 
corresponds to fi = 0.87. Neglecting the influence of the stationary contribution 
thus enhances the deviations from fj, = 1. 

The very good scaling of the 5 curves in Figs. 3.b, 3.c, involves 3 free pa- 
rameters: here A = 0.1, a = 0.02, /i = 0.97 or A = 0.1, a = 0.02, a = 2.2. They 
correspond to a full account of the results, with good coherence between TRM 
and x" experiments. Let us summarize our points: 

— The out-of-phase susceptibility x"(a;,t) is well represented by the sum of 
a stationary contribution, which varies like lj" (or (Ino;)^"""^''), plus a non- 
stationary contribution which varies as a negative power a: — 1 oiu.t (possibly 
with slowly varying corrections suggested by (26) above, but which are hardly 
visible on the x" data; see the rescaled x" in the inset of Fig.2). 



— The TRM decay curves are the sum of a stationary contribution oc t~°', 
plus an aging contribution. In the time range (f <C t«,) which is common to 
X" -measurements, this aging part can be approximated (in agreement with 
X") by a t/tyj-scalmg. The long time regime t >tyj seems however to prefer 
a sub-aging scaling of the type given in (17). 

The origin of this weak {1 — fi = 0.03) but insistant sub-aging behavior 
is an interesting point, still uncompletely understood. As mentionned above, 
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Fig. 3. d. Pull measured value of the TRM (as in Fig. lab), but vs the scaling variable 
A/To" defined in (29). 



from a theoretical point of view this would mean that the tq ^ limit of the 
underlying model does not exist. From a more physical point of view, several 
scenarios which might explain this effect can be considered. The simplest one 
concerns the effect of the field amplitude, which indeed (for larger fields) is 
known to suppress progressively the tw-eSect *(see [26] for details). However, 
the parameter /z seems to stick to a plateau value (< 1) for the explored low- 
field range of 10 — O.lOe (systematic scaling analysis with additive (rather than 
multiplicative) stationary corrections are however needed to confirm this point) . 

This sub-aging behavior might also be interpreted as a sign that aging is 
actually "interrupted" beyond very long, but finite, times [52]. For example, if 
the trapping time distribution is cut-off beyond a certain ergodic time terg (which 
itself depends on the subsystem), then for part of the subsystems ergodicity will 
be realized within the time of the measurement: their dynamics will no more 
depend on (interrupted aging). As shown in [52], this produces an effective 
sub-aging scaling very close to the /j. < 1 effect described here. The result is a 
value of a typical time terg, which must be understood as a crossover time scale, 
beyond which fi will further decrease to 0. In the analysis of [52], however, we 
had not properly taken into account the stationary contribution, which we have 
shown here to bring /x much closer to 1 (Figs. 3.b and 3.d). We had found tgrg 
of the order of 10^~^sec, which might therefore be underestimated. 



In a similar way, n is seen to decrease with the stress amplitude in polymer glasses 
[23] 
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5 Towards a Hierarchical Description of the Space of the 
Metastable States 

5.1 Temperature Variation Experiments 

Main Qualitative Features: ac Measurements. Until now, we have only 

presented resiilts which are obtained from aging experiments at constant temper- 
ature after the quench into the spin-glass phase. In another class of experiments, 
the temperature is varied during aging. In terms of thermally activated processes 
in a mountaneous free-energy landscape, one may expect to explore in more de- 
tail the various scales of the free-energy barriers which are involved in the slow 
dynamics. Thermal activation should - at first sight - be able to speed up or 
slow down the aging evolution. The conclusions of these experiments have been 
surprisingly instructive. 

An early result was obtained in [54] ; measuring the aging relaxation of x" in 
a CuMn spin glass, the authors observed that any step increase or decrease in 
temperature was causing an instantaneous increase of x" j followed by a slow de- 
crease. We may call this effect "restart of aging" , since the relaxation is renewed 
by the temperature change, which hereby produces a similar effect as obtained 
after the quench. This phenomenon reveals that the slow aging evolution towards 
equilibrium is significantly disturbed by relatively small temperature changes. 
Such a "chaotic dependence" of the equilibrium states on temperature has been 
predicted to occur as a consequence of frustration in [56] , where it is argued that 
the relatively small free-energy of an overturned region of spins (droplet) results 
from large cancellations at the surface of the droplet. These cancellations shoiild 
be very sensitive to temperature, hence the strong effect of a small temperature 
variation. 

When studied in more detail, the temperature variation experiments do not 
simply show a restart of aging upon any temperature change. We present in Fig. 
4 an experiment performed in such a way that the reaction of the x" -relaxation 
to either a decrease or an increase in temperature is very different [16]. The 
sample is first quenched from above Tg = 16.7 K to 12 K; due to aging, x" 
slowly relaxes. After 350 min, the temperature is decreased to 10 K. Despite the 
reduced thermal energy, the relaxation does not slow down, but restarts abruptly 
from a higher value, in agreement with [54]; this is the surprising chaotic-like 
effect, which looks as if the system was (at least partially) restarting aging from 
the quench. But, when after another 350 min at 10 K the sample is heated back 
to 12 K, the result is very different: x" resumes its slow relaxation from the value 
which had been reached before the temperature variation (see the inset of Fig. 
4), in a memory-like effect. 

There is no contradiction with the results in [54] . The experimental conditions 
of Fig. 4 are chosen here as the most illustrative of this twofold effect. Indeed, 
for smaller temperature variations, more intricate situations can be found [55]; 
aging at the lower temperature may partly contribute to aging at the higher 
temperature, even with a temporary incoherent transient, as emphasized e.g. in 
[17]. What we want to stress here is that a restart of aging is always observed 
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Fig. 4. Out of phase susceptibility x"('^i*a) of the CdCr1.7Ino.3S4 sample 
(Tg = 16.7K) during a temperature cycle. The frequency w is 0.01 Hz, and ta is 

the time elapsed from the quench. The inset shows that, despite the strong relaxation 
at 10 K, both parts at 12 K are in continuation of each other. 



after a temperature decrease; in contrast, what is found upon a temperature 
increase is a memory of previous aging at this higher temperature. On the one 
hand, if a long time has been spent previously at the higher temperature, as is the 
case in Fig. 4, only a weak relaxation is found; on the other hand, if the system 
has been directly quenched to a given temperature, heating up afterwards to a 
higher temperature - for which no memory of previous aging exists - will produce 
a strong relaxation at this temperature (see examples of various situations in 
[55]). 



Same Effects in TRM Experiments. The effect is quantitatively confirmed 
by measurements of the TRM- relaxation [25, 16]. The comparison with re- 
quires some care; the x" -relaxation directly shows on-line aging as a function of 
time, whereas in the TRM the effect of aging during tw is considered afterwards, 
during the relaxation which follows the field cut-off at tyj . The TRM curve shows 
the relaxation processes in a wide time window, and thus yields more extensive 
information than x" at a given frequency lu (which mainly reflects the processes 
of characteristic time ~ In Fig. 5. a, a negative temperature cycling has 

been applied during the waiting time. 

Prom the above x" results, one expects that aging processes be restarted 
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Fig. 5. a. Effect on the TRM relaxation of a negative temperature cycle 

(CdCri. 7/710. 3S4 sample, Tg = 16.7K). After waiting tw+ = 15mm at To 
(= 12K = 0.7Tg), the sample is cooled to To - AT for = lOOOmin, and then 
is heated back to To; after another tw+ = 15mm, the field is cut and the relaxation 
measured (at To). These relaxations (thin lines) are compared with normal ones, mear 
sured after waiting tw = 10,30, lOOorlOOOmm at constant To (bold lines). 



when going to Tq — AT; but for sufficiently large AT this evolution will be erased 

when coming back to Tq. This is what can be checked in Fig. 5. a for AT = IK: 
the procedure yields a relaxation curve which is exactly superimposed onto that 
obtained after simply waiting ty,+ + tyj+ = 30mm at constant Tq, aging during 
tw- ~ lOOOmin at 11 K has not contributed. Note that the equivalent normal 
curve has = 30mw = 2f^_|_, not 15mm (see Fig. 5a); the memory of the first 
aging stage has indeed been preserved. 

Intermediate AT values produce intermediate situations; for AT = 0.3K, 
the resulting curve is the same as that obtained after waiting tgff = lOOmin 
at constant Tq. Prom this example, we can work out a quantitative discussion 
of the effect. If we consider that the same characteristic free-energy barrier has 
been crossed by thermal activation during 

i) ty,+ + teff + tw+ = lOOmin at Tq and 

ii) {ty,+ at To) + (i^_ at To - AT) + (t^+ at To), 
then we can write 

(To-ZiT)ln^ =Toln^ , (31) 

To To 

which yields for the attempt time tq the unpleasant value of ^ 10~^^s. Obviously, 
the memory effect cannot be explained by the only thermal slowing down of 
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jumping processes over constant height barriers. 

In addition, simple thermal effects cannot be expected to explain the restart 
of aging, which is again evidenced in the positive cycling procedure of Fig. 5.b 
(all temperatures remain below Tg). 



Fig. 5. b. Effect on the TRM relaxation of a positive temperature cycle 
{CdCr1.-rIno.3S4 sample, Tg = 16.7K). A short heating cycle is applied after 970 min of 
waiting time at To; then one still waits 30 min before cutting the field and measuring 
the relaxation at To (thin lines). These relaxations are compared with normal ones, 
measured after waiting tw = 30 or 1000mm at constant To (bold lines). 

What appears in this procedure is the restart of aging due to the temperature 
decrease at the end of the cycle; for a sufficient AT (= 2.5K), the restart is so 
strong that the 970 min of previous aging are completely erased, as proved by 
the superposition of the AT = 2.5K curve with a normal = 30mm one. 

Again, for intermediate AT values, the effect is only partial; in its short-time 
part, the AT = IK curve sticks to the young ty, = 30mm-curve, whereas in 
its long-time part it goes closer to older ones. In contrast with the result in 
Fig. 5. a (negative T-cycling), the present procedure (positive T-cycling) yields 
curves which are not equivalent to a given waiting time at constant temperature. 
Clearly, the reason for that is that the positive temperature cycle ends by cooling 
down to To (which produces a partial restart of aging at To), whereas the negative 
temperature cycle ends by heating back to Tq (which lets the system retrieve 
the memory of the previous aging at Tq). 



AT (K) 




0.50 



t (min) 
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Thus, the same features of aging are observed in x" and TRM experiments 
[16, 25]. A chaotic nature [56] of the spin-glass phase appears when the tem- 
perature is decreased, the restart of aging processes being very similar to what 
initially happens after the quench. On the other hand, a memory effect is found 
when the temperature is raised back, and this memory effect goes far beyond 
what can be expected from thermal slowing down. 

A Hierarchical sketch. These effects have been interpreted in terms of a 
hierarchical organization of the metastable states as a function of temperature 
[25, 16]. The empirical picture is sketched in Fig. 6. 



Fig. 6. Schematic picture of the hierarchical structure of the metastable states as a 
function of temperature. 

During aging at a given T, the system samples the free-energy valleys (metasta- 
bles states) at a given level of a hierarchical tree; the restart of aging upon 
lowering the temperature is figured out as a subdvision of the free-energy val- 
leys into others, pictured as the nodes of a lower level of the tree, which thus 
develops multifurcating branches as the temperature decreases. The system is 
partially quenched since it must now search for equilibrium in a new, unexplored 
landscape, therefore aging (at least partially) restarts. The hierarchical picture 





1= 1, 2, 3 
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naturally provides us with the observed memory effect; when the temperature 
is raised back, the newly born valleys and barriers merge back to the previous 
free-energy landscape. Thus, aging at T — AT may have not contributed to the 
evolution at T, as far as all sub valleys which could be explored at T — AT origi- 
nate from unique valleys of the landscape at T (which is realized for large enough 
AT). 

The interpretation of these effects in a picture where aging is seen as the 
growth of compact, independent domains (droplets [6], domains [7]) remains, in 
our opinion, difficult ® . The droplets may be broken into smaller ones when the 
temperature is decreased, thus producing a restart of aging; but the memory 
effect implies that some information is kept somewhere about the stage of ag- 
ing which had been reached before decreasing the temperature. Thus, large-scale 
correlations should be kept untouched, while on the other hand they should 
apparently be destroyed. A way to satisfy these requirements could be to con- 
sider that the droplets should be fractal (non space filling) [8, 51], with a non 
trivial internal structure, in such a way that large droplets may contain smaller 
ones which can be activated independently. These droplets inside droplets might 
lead, under some energetic conditions, to a space-transcription of the hierarchical 
organization of the states. 

More Quantitatively: Rapid Growth of Free-Energy Barriers. The ob- 
served effect on aging of various thermal histories is in contradiction with thermal 

activation over constant height barriers, as evidenced above in (31). The picture 
can be made more quantitative along this same line (Fig. 5. a and 31), which 
we now further develop. In one experiment, the spin glass is aged during t~ at 
T— AT, and brought back to T before cutting the field and measuring the relax- 
ation. In another experiment, the spin glass is simply aged during t° at T, and 
the relaxation is measured at T. If both decay curves are the same in the whole 
time window, we may consider that the aging states reached in both procedures 
are the same, and that equivalent regions of both landscapes at T — AT and T 
have been explored. We can characterize each evolution by a typical height B of 
the maximum barrier which has been crossed in each case, namely 



Thanks to the identity of the TRM curves obtained in each procedure, one may 
consider that this quantifies the T- variation of the same barrier, which limits in 
both cases the same region of the phase space. The contradiction pointed out in 
(31) shows that B{T — AT) > B{T); in other words, the hierarchical picture of 
valleys bifurcating into valleys as the temperature decreases is supported by the 
observation of the growth of free-energy barriers. 

® See [25, 16]. However, some experiments have been analyzed along this line, with the 
introduction of long-time effects for the breakup of the domains [7, 17). 



B{T - AT) = (T - AT). In 



(32) 



B(r) = r.in^ . 

To 



(33) 
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An extensive series of TRM measurements on the AgMn sample, for multiple 
values oi t~, T and AT, has been performed [18]. In brief, the result is that 
the barriers are growing for decreasing temperatures below Tg, and in addition 
that their growth rate is so fast that at any T < Tg some of them should 
even diverge (see details in [18]). This provides us with an interesting link to the 
Parisi solution of the mean-field spin glass [4] . For decreasing temperatures, some 
barriers separating the mctastablc states arc diverging, transforming the valleys 
into pure states in the sense of the Parisi solution; the hierarchical structure 
of the valleys, deduced from the experiments, can thus be related to that of 
the pure states. Also, the picture which emerges from those results is that of a 
critical regime at any T < T^; in a sequence of micro-phase transitions starting 
at Tg, the spin- glass phase space continuously splits into nested and mutually 
inaccessible regions, within which non trivial dynamics takes place. 

5.2 Aging as a Random Walk: Traps on a Tree 

From a One-Level to a Multi-Level Tree. In its first stage [12], the trap 
model deals with all-connected traps, which can be called a "one-level tree" (see 
Sect. 3.3). The basic quantity which is calculated is the spin-spin correlation 

function CA{tw + t,ty^,), which can be related to the relaxation function if the 
fluctuation-dissipation theorem holds (in non-equilibrium, generalized forms of 
FDT might nevertheless hold [14], see Sect. 3.2). The aging part Ma of the 
TRM-rclaxation can thus be estimated from the decay of the aging correlation 
function, which is proportional (by a factor qea) to the probability n{t, tw) that 
the system has not jumped out of a f^-trap at -|- 1: 



The probability iJ contains all the information on possible jumps from trap to 
trap, and thus represents the aging dynamics; it should tend to 1 when t^ goes to 
infinity, for any finite t. In this limit, equilibrium dynamics is recovered, since no 
jump occurs; the proportionality factor qEA in (34) describes this "bottom of the 
traps" dynamics (c/. with (14) in Sect. 3. 2). From its definition, it represents the 
overlap between the various configurations which constitute the bottom of the 
traps; in a TRM experiment, it is this finite fraction of the initial magnetization 
which should be found in the ideal limit of infinite tw 

Two asymptotic behaviors of 77 (t, i^) (and thus of the TRM) are calculated 
in a "multi-level" version of the trap model [13]: 



The predicted shape of the aging TRM is thus a constant minus a power 
law of exponent 1 — Xm for the short times [t 4C tw), and on the other hand 



M^(i 4- tw,tw) ~ CA{tw + t, tw) = qEA-n(t, tw) 



(34) 




(35) 
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a simple power law of exponent xi at long times {t ty,). In the simple one- 
level case of all-connected traps with trapping times distributed with a given x 
(27), one has xi = xm = x, which does not yield a realistic TRM-decay shape. 
However, a satisfactory fit to the TRM measm'cments can be obtained with only 
two values of x [13, 52]; one finds from the initial part of the TRM (or from 
x" measurements) xm ~ 0.65 — 0.8 and from its late part xi = 0.05 — 0.35 (for 
several samples and different temperatures). In Fig. 7, we show the asymptotic 
behaviors corresponding to (35) for a typical TRM curve, 1—Xm being equivalent 
to the exponent in (22) if h{t) is a power law. 
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Fig. 7. Aging part of one of the above TRM curves (Fig. 1 and 3). Both asymptotic 
limits t <C t„: and t ^ have been fitted to (35). The variation of the effective index x 
along the different time regimes suggests a more complex than one-level tree structure. 



Thus, if one thinks in terms of a one- level picture with no further assump- 
tion, it appears from the shape of the TRM itself that the effective x which 
labels the trapping time distribution (27) is different in the short and long-time 
regimes. Namely, x is closer to 1 for i <C [x = xm), and decreases for t^^tw 
{x = x\). This simple observation has been translated in terms of a hierarchical 
organization of the traps (traps inside traps); the model of all-connected traps 
[12] (one-level tree) has thus been generalized in [13], of which we now extract 
the main points. 

The one-level tree is a distribution of traps of index say x = xi (< 1); 
following the construction of the Parisi tree of states [4], it can be extended 
into an n-level tree by multifurcating each trap into others of index x^ > xi, 
which themselves subdivide into others, until some final index a;„ > 1 at the 
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end-branches of the tree. At each level i, the Xj-distribution of trapping times 
corresponds [4, 12] to an Xj-dependent distribution of free energies, and this 
yields for the overlap between the corresponding states a value qi{xi) which is 
the inverse of the Parisi order parameter q(x) [4], thus (jj+i > (jj. Now, the states 
close to the end-branches (say the lower part of the tree) with a; > 1 can all be 
visited in finite times, since for a; > 1 the mean value of the distribution (27) is 
finite: the corresponding dynamics is stationary. 

Aging phenomena appear when x < 1, i.e. above a given level c such that 
Xc = 1. The faster processes of the aging regime, which are seen in the short-time 
part of tlic^ TRM, occur among the states which are close to each other in the 
hierarchical geometry (large overlap), that is which are related by a tree- node 
at a level very does to c; therefore they correspond to x close to 1, as is indeed 
suggested by the exponent of the power law behavior at the beginning of the 
TRM. As time elapses, more distant states (of smaller overlap) can be explored, 
and the corresponding transitions imply passing higher nodes in the tree, of 
smaller index x, in agreement with the smaller exponent of the power law in the 
long-time part of the TRM. 

The model can be solved for an arbitrary number of levels; but the simpler 
case of a two-level tree has been computed and fitted to the TRM experiments 
[13] . The fits are of a very good quality for both insulating and metallic samples: 
for a given t^, the shape of the decay curve is quite well reproduced over the 
whole t scale. If the number of metastable states at each level is infinite, however, 
the predicted influence of t^ is precisely a t/tw scaling (full aging), as in [12], 
and this does not satisfactorily account for the measured tw-dependence, which 
corresponds to a sub-aging behavior (see the discussion in Sect. 4). 

Effect of Temperature Changes. Now, the effect of the temperature on this 
tree-like organization of the metastable states can also be discussed. From the 
fit of the TRM-data [52], the XM-index (e.g. from the short-time part) shows 
a tendency to increase towards 1 as T approaches Tg-. this suggests that the 
spin-glass transition be closer to the REM scenario [50], where x = T/Tg, than 
to the mean-field equilibrium scenario [4] where the Parisi parameter x rather 
goes to zero at Tg. This leads us to a possible interpretation of the empirical 
hierarchical structure versus temperature in terms of the Parisi-like tree of traps 
[13]. 

Tg is characterized by a certain c- level of the tree where Xc = l: the tempera- 
ture dependence of x means that, when the temperature is lowered, this is now a 
lower (multifurcated) level of the tree which corresponds to x = 1 and limits the 
stationary dynamics. Thus, aging dynamics restarts among newly born states, as 
observed in the experiments and interpreted along the hierarchical scheme of Fig. 
6. Conversely, when the temperature is raised back, the lower aging level comes 
back to equilibrium, and aging continues at the upper level, resuming from its 
previous stage of evolution. A more detailed comparison between theory and ex- 
periment is however desirable. Note finally that a hierarchical decoupling of time 
scales has also been argued on the basis of "second" noise spectra by Weissman 
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and collaborators [58]. 

The influence of temperature is thus crucial in these pictures where activated 
processes play the dominant part in aging. Conversely, temperature is somewhat 
irrelevant for the aging dynamics of the mean-field models described above (at 
least those corresponding to single-scale dynamics). For example, the effect of 
temperature variations has been explicitely computed within the spherical p = 2 
model [57], with the result that nothing comparable to experiments happens. In 
this respect, aging in these single scale models is again very similar to simple 
domain growth in a ferromagnet [21], which is rather insensitive to temperature. 
It would be very interesting to understand how temperature changes affect the 
dynamics of a multiscale model, such as the SK model. 



6 Conclusions 

In this paper, we have presented a survey of the experimental results concerning 
aging phenomena in spin glasses, focusing on the description of magnetization 
relaxation and low- frequency out-of-phase susceptibility measurements. We have 
discussed in some details how far two theoretical approaches of out-of-equilibrium 
effects, namely the trap model [12, 13, 52] and the microscopic approach [14, 15], 
can succeed in describing different aspects of the dynamics. 

The relaxation of the out-of-phase susceptibility x" towards a non-zero value 
indicates the presence of stationary dynamics. This same dynamics can be found, 
although less obviously, in the short-time part (compared to the waiting time tw) 
of the decay of the thermo-remanent magnetization (TRM). In fact, an additive 
combination of a stationary and an aging regime in the auto-correlation (and 
thus also in the TRM and in x") follows from the solution of some mean-field 
models [14, 15, 38, 39, 57]. This striking similarity of mean-field and real spin 
glasses raises the question of the nature of the slow dynamics in these systems. 
While one is used to think of thermal activation in mountaneous landscapes as 
the source of slow processes, the mean-field models now provide us with aging 
phenomena which are due to the flatness of large regions in the phase space [21], 
with no crucial role played by the temperature. 

For the sake of a coherent description of both TRM and x" data, one should 
therefore extract the same stationary contribution from all results, which we 
have made here in an additive way. The remaining aging contribution presents, 
for the TRM, systematic departures from a "full aging" situation of a pure t/t^ 
scaling ("sub-aging") [3, 28]. When the stationary dynamics of the TRM is ad- 
ditively accounted for, the departures from full aging become less pronounced, 
but still remain. For x" , whose aging regime only overlaps that of the TRM on a 
limited range, both full aging and sub-aging scalings remain compatible with the 
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data. Using some recent analytical developments of the microscopic approach of 
spin-glass dynamics [15], we have shown that a sub-aging behavior will appear 
under some general conditions. Thus, the microscopic theory can now account 
for the scaling functions which had been postulated in the past by the experi- 
mentalists on phenomenological grounds [3, 28]. 

The experiments have shown that small temperature variations have a strong 
effect on aging phenomena [16, 17, 18]. A temperature decrease restarts the evo- 
lution, whereas raising back the temperature lets the system retrieve its previous 
stage of aging at this same temperature. These results have been interpreted in 
terms of a hierarchical organization of the mctastable states as a function of 
temperature, in which the valleys of the free-energy landscape subdivide into 
others for decreasing temperatures [16, 18]. The mean- field models have not yet 
brought conclusive results on this question. The trap model [12], which provides 
us with a stochastic picture of aging, has recently been extended [13] in a way 
which sheds some light on the temperature variation experiments. 

The simple picture of a random walk among all-connected traps (one-level 
tree of states) did not take into account the existence of stationary dynamics, 
which in this language is a "bottom-of-the-traps" dynamics. For this reason, 
traps must have an internal structure. On the other hand, the TRM shape is 
related to the index x of the trapping time distribution. The comparison of the 
experimental shapes with a one-level tree model shows that the "effective x" 
systematically varies along the curve; it is closer to 1 in the t <C tu.. regime, and 
smaller in the t ^ tyj regime [12, 52, 13]. Thus, the traps which are explored in 
the various time regimes are not connected in the same way. 

These results are fairly well understood within a picture of hierarchically 
connected traps, in a multi-level tree geometry where x varies from one level 
to the other. The limit value x = 1 sets the tree level which is at the border of 
stationary and non-stationary dynamics. The tendency of x to increase when the 
temperature approaches Tg [52] indicates that this border x = 1 level changes 
with temperature; the resulting model of a tree-like organization where the limit 
of equilibrium dynamics changes with temperature [13] is now very close to the 
empirical picture of a hierarchy of metastable states as a function of temperature 
which had been drawn from the experiments [16]. 

This recent development, together with the predictions of slow dynamics and 
aging from microscopic models, have created a very stimulating atmosphere. 
They let us expect that an even tighter interaction among theoreticians and 
experimentalists will bring us closer to a more complete understanding of the 
physics of disordered systems. 
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